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Abstract 

The two-state model of stochastic resonance is extended to a chain of coupled 

two-state elements governed by the dynamics of Glauber's stochastic Ising 

model. Appropriate assumptions on the model parameters turn the chain into 

a prototype system of coupled stochastic resonators. In a weak-signal limit 

analytical expressions are derived for the spectral power amplification and 

the signal-to-noise ratio of a two-state element embedded into the chain. The 

effect of the coupling between the elements on both quantities is analysed and 

array-enhanced stochastic resonance is established for pure as well as noisy 

periodic signals. The coupling-induced improvement of the SNR compared to 

an uncoupled element is shown to be limited by a factor four which is only 

reached for vanishing input noise. 
PACS numbers: 5.40. +j, 05.50+q 
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I. INTRODUCTION 



The essential point of stochastic resonance [1] is the following: If a stochastic resonator is 
subjected to an external influence, commonly referred to as the signal, its response exhibits 
signal features, which are most pronounced at a certain level of noise present in the system. 
The fascinating aspect of stochastic resonance is therefore that, counter- intuitively, an in- 
creasing noise level does not steadily deteriorate the transmission of the signal through the 
resonator. On the contrary, noise may rather be used to optimise the transmission process 
[2]- 

The implications of these effects are intensively investigated. In biological systems, for 
example, stochastic resonance apparently plays a role in the neural transmission of informa- 
tion [3] . From a technical point of view one might possibly see the emergence of a novel type 
of detector which incorporates an optimal amount of noise to perform best. In this regard 
superconducting quantum interference devices (SQUIDs) have been studied [4]. They may 
be made to detect weak magnetic fields reasonably well without the usual costly shielding 
from environmental noise. 

In recent years it was shown that the performance of a single stochastic resonator can 
be enhanced, if it is embedded into an ensemble of other stochastic resonators which are 
properly coupled [5] [6] [7] [8] [9]: Compared to being operated isolated the response of 
the resonator to the signal increases within the coupled ensemble. However, if the coupling 
becomes too strong this response was found to deteriorate again. Thus in an ensemble of 
stochastic resonators the coupling strength turns out to be a second design parameter: Apart 
from the noise level it can be tuned to achieve an optimal performance of the embedded 
resonator. This effect was called array- enhanced stochastic resonance [8]. It will possibly 
find technical exploitation and might also be relevant to biological systems, for example, 
coupled neurons. 

The aim of the present paper is a further investigation of the phenomena related to 
array-enhanced stochastic resonance. To this end we study a simple prototype system of 
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coupled two-state stochastic resonators under periodic modulation. The model we propose 
allows to analytically calculate the weak-signal limits of two prominent stochastic-resonance 
quantities, spectral power amplification (SPA) and signal-to-noise ratio (SNR), respectively 
Following the concept of array-enhanced stochastic resonance both quantities refer to the 
response of a resonator as part of the coupled ensemble. The impact of the coupling on both 
response measures will be studied. 

There is a close link of the present analysis to the work presented in [5] and [8], which 
to our knowledge are the most comprehensive in the field. However, compared to [5], where 
the SPA of a system of globally coupled bistable elements was studied analytically, too, 
the present approach offers a more refined description, since it allows to avoid a mean-field 
approximation. Compared to [8], on the other hand, where the SNR of a chain of locally 
coupled elements was investigated in a (rather extensive) simulation, the SNR can now be 
obtained analytically. For a quick review of the results the reader is referred to the 3rd and 
4th passage of the summary section. 

From a conceptual point of view the model proposed here can be seen as an extension 
of the two-state model developed by McNamara and Wiesenfeld [10] to study stochastic 
resonance in noisy bistable systems. Instead of considering individual two-state elements 
as in [10] these elements are now arranged in a chain. A simple next-neighbour interaction 
is introduced which brings the model close to the system studied in [8]. The interaction is 
chosen in a way that the resulting evolution of the elements is given by Glauber's stochastic 
Ising model [11]. A detailed description of this approach is given in the next section which 
also provides the necessary background. We note that the connection to the Glauber model is 
made for mathematical convenience. We are not concerned with the observation of stochastic 
resonance in Ising systems which was the central theme in [12] and [13]. Our purpose will 
also lead to the assumption of Arrhenius-type transition rates of uncoupled unmodulated 
elements which is unusual in the context of Ising systems. First results of the present analysis 
had been published in [14]. 

A further interesting approach to arrays of stochastic resonators is c entered around a 
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response of a more collective nature. Here stochastic resonance is studied in the summed 
output of N resonators which need not necessarily be coupled [15] [16] [17]. Recently, it was 
shown analytically that the SNR of this output approaches the input SNR, if a sufficiently 
large number of uncoupled resonators is used [18]. Whether in some way coupling may 
still be beneficial for the collective response of the present model will be the subject of 
a future study. Only the summed output of a very large number of resonators shall be 
briefly discussed here. It will be shown that in this case the SNR always deteriorates under 
coupling. This indicates that the coupling-induced improvement of the performance of 
stochastic resonators associated with array-enhanced stochastic resonance is a local rather 
than a global effect. 

II. FROM SINGLE TO COUPLED TWO-STATE STOCHASTIC RESONATORS 

An early theory of stochastic resonance in noisy bistable systems with one (generally 
continuous) variable was worked out by McNamara and Wiesenfeld [10]. They studied the 
effect of a periodic modulation of these systems in a way which is independent of the precise 
dynamics involved. To this end the behaviour of the bistable system was approximated 
by a random telegraph process. For simplicity this process was taken to be symmetric, 
randomly switching between two states a = ±c. For convenience we shall consider here 
c = 1. The probabilities p(cr) to find the process in state a satisfy p(cr) + p(—cr) = 1 and 
their time-evolution is governed by 



where W(a) denotes the rate of the transition a — > —a. These rates must be extracted from 
the precise dynamics at hand. To build a theory of stochastic resonance they were assumed 
to be approximately given by [10] 



p(a) 



.p(- a ) = w(-a)p(-a) - W(a)p(a), 



(1) 




(2) 
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with 5 being the small parameter of the theory. Obviously, at this level the periodic modu- 
lation of the bistable system with (angular) frequency cu and phase <f> is taken into account 
by the cosine term, only. 

This simple model allows to calculate the long-time limit of the power spectrum of the 
random telegraph process averaged over a uniformly distributed phase 0. The spectrum 
consists of a continuous part and a delta function at modulation frequency uj. The contin- 
uous part of the spectrum is the Lorentzian of the unperturbed process times a frequency- 
dependent prefactor. The latter is close to 1 and governs the modulation-induced transfer 
of broad-band power to the delta peak in the spectrum. 

From this spectrum McNamara and Wiesenfeld obtained their central theoretical result 
with respect to stochastic resonance: An exact analytical expression for the signal-to-noise 
ratio (SNR) of the two-state process. This SNR is defined as the ratio of the weight of the 
delta function to the continuous part of the spectrum at modulation frequency. In view of 
the present work we shall neglect the signal-induced suppression of the continuous part of 
the spectrum and consider the linear response approximation of the McNamara- Wiesenfeld 
result with respect to the small parameter 5. Doing so one finds for the SNR 



If the approximation (2) is performed on a particular noisy bistable system subject to 
a periodic modulation, the parameters a and 5 become functions of noise strength and 
modulation amplitude, respectively. The dependence of the SNR on the noise intensity 
can then be studied and the occurrence of stochastic resonance may be established for that 
particular system. 

As an example the overdamped double- well system was given in [10]. The model equation 
reads 



where ^(t) is Gaussian white noise with (£(£)) = and (£(t)£(t+r)) = 5(t). Using a modified 
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(4) 



Kramers formula for the transition rates valid for sufficiently low modulation frequencies 
McNamara and Wiesenfeld found 



The central idea of the present work is to pass from the single bistable element to a set of 
coupled bistable elements by extending the two-state model of McNamara and Wiesenfeld. 
To this end we consider a chain of two-state elements which for convenience is taken to be of 
infinite length. If the elements interact the simple gain- loss balance (1) has to be modified 
to describe the evolution of the set of probabilities p(a, t) to find the chain in a particular 

configuration a = ( , <Jk-i, ^k, <?k+i, ) at time t. Introducing a formal operator F k 

defined for any function f(cr k ) by F k f(a k ) = /(— a k ) the new gain-loss balance reads 



With the initial condition p(a, t ) = 5 5 ^ the p(a, t) are interpreted as the transition proba- 
bilities p(a,t\a ,t ) of a Markovian process with infinitely many discrete components. The 
statistics of the process is fully determined by (6). 

To introduce interactions we assume that the transition rates (2) of an element depend on 
the states of its next neighbours. A simple choice for this coupling controlled by a parameter 
7 is 



where W(ai) is the McNamara- Wiesenfeld rate (2). With positive 7 neighbouring elements 
prefer to be in the same state, whereas they tend to be in opposite states, if 7 is negative. 
Both tendencies grow with growing coupling strength |7|. To avoid negative transition 
rates I7I < 1 has to be required. We note that the assumed type of coupling is thus able 
to model ferromagnetic-type interactions for which a coupling-induced improvement of the 
performance of stochastic resonators was found in [5] and [8]. For positive 7 the present 




(5) 



p(a) = £(F fe - l)W k (a k )p(a). 



(6) 



k 



W^) = W(<Ti) 1 - ^(<7i_i + <T i+1 )<Ti 



(7) 
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model is particularly close to the model simulated in [8], since both are linear arrays of next- 
neighbour-coupled elements. For convenience we shall often use the terms ferromagnetic or 
antiferromagnetic coupling instead of coupling with positive or negative 7. 

Of course, other types of transition rates (7) could be considered as well. The advantage 
of the present choice is that most of the relevant stochastic properties of the resulting model 
are already known. They were studied by Glauber who introduced this model as a stochastic 
form of the Ising model [11]. (We note that in the Glauber model the term <5cos(u;t + 0) in 
(7) is replaced by a general time-dependent parameter (3.) As in the McNamara-Wiesenfeld 
model it is again possible to find an analytical expression for the SNR in leading order of the 
modulation parameter 5 which is the subject of the next section. The respective analysis is 
essentially an exploitation of Glauber's work. 

It may be possible to reduce the dynamics of some coupled systems to the two-state 
model given by (6) and (7) and hence to directly employ the SNR formula derived in the 
next section. However, this approach is certainly more involved than the rate expansion (2) 
needed to make use of the McNamara-Wiesenfeld SNR. We shall not pursue this more general 
aspect of the model in further detail. Instead we are interested in devising a simple prototype 
system which allows to study the impact of the coupling on the stochastic-resonance effect. 
To this end we are looking for simple assumptions on the dependence of the model parameter 
a, 7, and 5 on some noise intensity and signal amplitude, respectively. We proceed in two 
steps: 

First, we retain the Glauber's original relation of the present model to the Ising model 
given by the Hamiltonian 



This relation is defined by an assumption on the probability p(<7 k ) to find the k-th element 
in state a k , if all other elements are fixed: In an adiabatic limit the p(<Jk) of both models 
are to be identical. This implies that 



(8) 



k k 



P(-Cfc) 



W k (-a k ) 
W k (a k ) 



= exp 
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holds which allows to establish a relation between the parameters of both models. One finds 
7 = tanh(2 J/kT) and <5 cos(W + 0) = tanh(/x# j kT) , respectively [11]. The role of the noise 
intensity is thus played by the temperature T. To skip unnecessary parameters we set k — 1 
and \i = 1. In the present paper we assume that H, which shall be referred to as the signal, 
is given by 

H = H cos(ut + (f)). (10) 

For weak amplitudes H <T this results in S — H /T. (In section V we shall also consider 
periodic signals (10) with an additional noisy background.) 

Second, a plausible assumption has to be made on the temperature dependence of the 
time-scale parameter a which is not effected by requiring (9). It is assumed that it retains the 
qualitative dependence (5) found by McNamara and Wiesenfeld for the double-well system, 
i.e., a = tto ex p( _ 1/^) with properly scaled temperature T. This dependence is typical for 
many rate processes [19]. An appropriate temperature dependence of a is also necessary to 
turn the elements into stochastic resonators for vanishing coupling. For simplicity a = 1 is 
considered. 

The prototype system we are looking for is thus specified by (6) and (7) together with 

a = exp (-^) , 7 = tanh , 6=^-. (11) 

We emphasise that this system is not an approximation of the Ising system (8) in the same 
way as the double-well system (4) is approximated by (1),(2), and (5), respectively: The 
Ising system and the stochastic Glauber model are based on totally different dynamical 
concepts. Nevertheless the devised prototype system appears to capture in a simple way the 
essence of the effect of coupling and forcing on coupled continuous-variable systems like the 
double well. 

To motivate this, let AU(a) denote the height of the barrier to be surmounted in order 
to escape from state a. In a very simple approximation the effect of next-neighbour coupling 
and forcing on the barrier height could be modelled as 
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AU(a k ) = AU + AU aM (a k ), 



(12) 



AU add (a k ) = Ja k (a k+1 + <7 fe _i) + H a k cos(wt + 0). 



Comparing this to (8) one finds that within this simple picture the effect of coupling and 
forcing on the barrier heights is the same as their effect on the energy levels of the Ising 
model. (We note that for coupled double-well systems the chosen approximation would only 
hold for weak coupling. For strong coupling the systems may no longer be bistable and the 
entire concept of a simple-minded two-state approximation breaks down.) 

Within an adiabatic limit and with some noise strength T the transition rate W(a k ) of 
the barrier system would roughly read 



From here rates of the form (7) together with parameters (11) can be obtained in two steps. 
First, the factor related to AU is kept as the parameter a with a respective choice of AU . 
Second, the factor related to AU a dd, which describes the effect of forcing and coupling, is 
replaced by the respective term of (7), which is easier to handle. A connection between both 
types of rates is again made by requiring that they lead to identical stationary distributions 
in the sense of (9). 

A particular close relation between the Glauber-dynamics model and the barrier system 
should thus arise within the regime of quasistatic response. Nevertheless we shall not restrict 
the following investigation of the simple prototype system to the quasistatic regime. Whether 
or not the suggested model provides insights which are also relevant to more realistic settings 
and possibly helpful in the design of coupled stochastic resonators remains to be seen. 




(13) 



III. ANALYTICAL SOLUTIONS 



The aim of this section is to calculate the SNR of the response of a two-state element 
embedded into the chain of coupled elements given by (6) and (7). This calculation will 

9 



be done in terms of the parameters a, 7, and 5 and is essentially based on Glauber's work. 
We shall only take account of terms in leading order of the modulation amplitude 5. The 
resulting SNR will then be an extension of the simplified McNamara-Wiesenfeld result (3). 
At the end of the section the general result will be applied to the prototype system specified 
by (11). In addition to the SNR we shall also give an analytical expression for the so-called 
spectral power amplification (SPA). 

Following Glauber we start by calculating the averages (ak(t)). From (6) one derives 
d(ak)/dt = —2(a k Wk(<Jk))- It results in 



with Tij(t) = (<Ji(t)<Tj(t)). A closed set of equations for the (<Jk(t)) is found by linearising 
(14) in 5. In this case the r iy j can be taken from the unperturbed model (5 = 0). After long 



The resulting long-time-limit set of equations can be further simplified. Since all elements 
are forced uniformly, a particular element cannot be distinguished from any other, once the 
initial distribution has been forgotten. Hence, all elements will have identical statistics and 
the indices in (14) can be skipped. Together with the previous assumption one obtains 



Obviously, the long-time dynamics of the average (a) is identical to the long-time dy- 
namics of the average of an uncoupled element with rescaled relaxation rate a(l — 7) and 
rescaled modulation term a\J\ — 7 2 5. The resulting long-time limit of the averaged state 
simply reads 



d(<7k) 
d(at) 



-fa) + 2 K^fc-i) + fa+i)] 

+ 1 - ^(^fe-i.fe + r k ,k+i) 5 cos(ut + <fr) 



(14) 



times these r i: j read = rj l ^ with i] = 7 1 (1 — ^1 — 7 2 ) [11]. 




(15) 



(<r(t)) = q cos (cut + (f) + , 



(16) 
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where 



q s is the response to static signals {uj = 0): 



= 5 



1+7 



(17) 



1-7 



Formally this static response might even exceed q s — 1. We conclude that in this case 5 is 
not sufficiently small and the linearisation of (14) is no longer justified. Later on it will be 
shown that this does not amount to a considerable restriction. 

The power spectrum of an element is determined from the average {a k {t)a k {t + r)). Like 
in [10] its t-dependence is removed by averaging (.)</, over an uniformly distributed initial 
phase (j) °f the modulation term in (7). Introducing the correlation function c kk (t,r) = 
([a k (t) - (a k (t))][a k (t + r) - (a k (t + r))]) one finds 



The second term on the r.h.s. contributes a delta function with weight nq 2 at signal frequency 
uj to the (one-sided) power spectrum as defined below. The first term, on the other hand, 
forms the continuous part of the spectrum. 

For the purpose of calculating the SNR in leading order of the modulation amplitude 5 it 
is sufficient to approximate this continuous part by the power spectrum of the unperturbed 
model. This implies that we neglect any possible signal-induced transformation of the con- 
tinuous part of the spectrum as we did to obtain the simplified McNamara-Wiesenfeld result 
(3). For the unperturbed model the long-time correlation function c(r) is not effected by the 
average (.)^. It also no longer depends on the index of the element. According to Glauber 
c(t) reads 



(Mt)crk(t + t)))^ = (c kk (t, r)) + — cos(cjt). 



(18) 



c (r) = e -«M £ n^Ua^Tl) 



(19) 



n=— oo 



with modified Bessel functions 7 n (a7|r|). Using the relation [20] 




(20) 
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which holds for Re(v) > — 1 and Re(a) > \Re(b)\ the power spectrum s(Q) = 
2 J +o ° cos (fir) c(r)dT is calculated. For the one-sided spectrum defined by S(Q) = s(Q) + 
s(— Q) at Q > we eventually obtain 

with si = \JJa + ifi) 2 — («7) 2 and S2 = 0:7(0; + «f2 + si) -1 . 

Now the SNR of the response of a two-state element embedded into the chain can be 
calculated, which is the central result of this section. It reads 

R* = ^fr (22) 
S{uj) 

For vanishing coupling it reduces to the simplified McNamara-Wiesenfeld result (3). 

In addition to the effect of the coupling on the SNR we also wish to study its impact on 
the delta peak in the power spectrum. In this regard the spectral power amplification (SPA) 
is a convenient measure (cf. e.g. Jung in [2]). It is defined as the ratio of power contained in 
the signal peaks of output to input spectrum, i.e., the ratio of the weights of the respective 
delta functions. Because of its dependence on the input spectrum the SPA depends on the 
precise dynamics to be modelled by the two-state chain. 

For the prototype system specified by (11) the SPA reads p = (q/H ) 2 , or explicitly 



/ 

P = Ps 

V 



1 + 7 7^2 , (23) 



(l-tanh(f))' 



Ps = y2^(y)' (24) 

where p s is the SPA of static signals. Obviously, the SPA does no longer depend on the 
signal amplitude H . For convenience we also remove the ifo-dependence of the SNR by 
considering 

K ~ H$S(uj) ~ S(u>Y [ ] 
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For the prototype system the parameter 77 involved in the unperturbed spectrum (21) sim- 
plifies to 7] = tanh( J/T) [11]. 

Both quantities, SPA (23) and rescaled SNR (25), only depend on three parameters: 
temperature T, coupling strength J, and signal frequency u>. Restrictions arise from H <C T, 
which lead to a simple expression for 5, as well as from the linearisation of (14), which results 
in an upper bound on 5 as discussed in connection with (17). Together one finds 



It implies that finite Hq place a lower bound on T and an upper bound on J. However, these 
restrictions are rather weak: H can be made arbitrarily small because its size is immaterial 
within the present weak-signal approximation. 

We are now in a position to study the impact of the coupling on the response measures 
SPA and SNR of a single resonator embedded into the chain. This will be the subject of the 
following two sections. 

IV. COUPLING AND SPECTRAL POWER AMPLIFICATION 

The SPA (23) has an unique maximum over temperature T and coupling parameter J for 
any time- dependent signal (10). The maximum SPA is obtained for a frequency-dependent 
value J max which is always positive. In other words: A properly tuned ferromagnetic- 
type coupling yields the best SPA performance of the two-state resonator element which is 
embedded into the chain. 

For our simple model this maximum can be studied analytically. The partial derivatives 
of the SPA with respect to T and J are found to vanish at pairs (T max , J max ) given by 



H Q < min(T,Texp(-2J/T)). 



(26) 



max 



(T max /4)ln(2T r 



max 



1) 



(27) 



and 




(28) 
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respectively. For J max to be real and finite, T max > 1/2 has to hold. With this restriction 
(28) has a unique solution 1/2 < T max < 1 for any given uo > 0. The corresponding J max 
is found via (27). Within the given range of T max one can easily check on (27) that J max is 
indeed always positive. An inspection of the second derivatives of the SPA finally reveals 
that there is a maximum at (T max , J max )- Its peak value is implicitly given by 

1 — T 

o - 1^ (29) 

1 max\ Z1 max 1 ) 

From (27), (28), and (29) the following features of the maximum can be derived which 
are shown in Fig. 1 : Tuning the signal frequency uo from very large to vanishing small values 
the temperature T max falls from 1 to 1/2. At the same time, the coupling strength J max as 
well as the peak height p max increase from vanishing small to very large values. In addition, 
an analysis of the curvatures pjj and Ptt, both expressed in terms of T max , shows that 
the sharpness of the maximum grows as its height increases. This indicates that to achieve 
optimal performance the system parameters T and J have to be tuned with increasing 
accuracy as the signal frequency decreases. 

Fig.l also includes a comparison to the SPA of an uncoupled element. At given uo > 
this SPA has a maximum located at a temperature T , implicitly given by 

2 / 2 \ T 



on the interval < T < 1, whereby its peak value is found to be p = (1 — T )/Tq. It can be 
shown from (28) and (30) that T max always exceeds T . The fact that J max was found to be 
never zero also implies that p max always exceeds p . However, looking at Fig.l it is obvious 
that at frequencies uo > 1 both compared quantities of the coupled element approach those 
of the uncoupled one and J ma x approaches zero. The coupling-induced increase in the SPA 
is thus vanishingly small at sufficiently high frequencies of the signal. 

This comparison allows to distinguish between two different types of SPA behaviour: At 
low frequencies the SPA is enhanced under ferromagnetic coupling whereas at high frequen- 
cies it is basicly not, although a tiny increase still occurs. Both situations are illustrated 
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at selected frequencies in Fig. 2 and Fig.3, respectively. The low frequency SPA exactly re- 
produces the qualitative effect found analytically by Jung et al. [5] in a system of globally 
interacting elements. The high frequency behaviour, on the other hand, where the SPA does 
practically not increase under coupling, has to our knowledge not been reported before. 

In both cases Fig.2 and Fig.3 clearly demonstrate that at fixed coupling parameter J 
stochastic resonance occurs: The SPA has a maximum over temperature T which, unfor- 
tunately, cannot be established analytically. Fig. 2 (bottom) shows that while ferromagnetic 
coupling improves the SPA at any fixed temperature T this improvement is lost, if the cou- 
pling becomes too strong. There is thus an optimal coupling strength J op t(T) for every given 
temperature, which can even be exactly calculated at any uj > 0. One finds 



which is always positive. Since the SPA does not have a further extremum over J, it implies 
that antiferromagnetic coupling always decreases the SPA. Moreover, J opt (T) increases with 
growing temperature T and decreases as the signal frequency grows. In Fig.2 the SPA at 
the optimal coupling strength is included (dash-dotted curves). In Fig. 3 (bottom), where the 
high-frequency SPA is illustrated, the coupling-induced improvement of the SPA is hardly 
detectable and J opt (T) is almost zero. 

Considering the SPA of time-independent signals (uj = 0) given by (24) one finds that 
it has a maximum over temperature for antiferromagnetic coupling J < 0, only, whereby it 
decreases as the coupling strength |J| grows. For J > the SPA increases with increasing 
J as well as with decreasing T. It formally diverges for T — > and J — > oo, respectively. In 
both cases the weak-signal limit (26) breaks down. 

In general, the SPA (23) is given by its static value (24) times a dynamical factor (l+rf 2 ) -1 
(cf.(16)). Here d is the ratio d = uj/(a(l — 7)) of signal frequency to long-time relaxation rate, 
which also governs the phase shift ip in (16). With growing d the elements gradually loose 
their ability to follow the signal: The SPA weakens and the phase shift grows. This effect 
occurs, for example, if the signal frequency uj increases. Subsequently, the SPA decreases 




(31) 
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with growing u. 

The impact of d is also responsible for the occurrence of stochastic resonance and optimal 
coupling in the SPA. Here a decrease of the long-time relaxation rate a(l — 7) plays the 
crucial role: The elements dynamics slows down as J increases or as T decreases. Thus 
this slow-down occurs whenever the static SPA grows. Hence, the dynamical factor always 
counteracts the static SPA as T or J is changed. Eventually, the increase of the monotonous 
static SPA is outperformed by the decrease of the dynamical factor which results into a 
maximum of the SPA over T and J, respectively. In other words, stochastic resonance and 
optimal coupling occur. 

As illustrated in Fig.2 the SPA can be seen as a transition between two limits, the static 
SPA given by p s = (q s /H ) 2 (d <C 1, dashed curves) and p s /d 2 (d 3> 1, dotted curves), 
respectively. Since both limits intersect at d — 1, i.e. at oj = a(l — 7), the plots nicely show 
the well-known approximate matching of time scales at the SPA peak. This matching does 
not only occur over T (Fig.2, top), but over J (Fig.2, bottom), too. If the signal frequency 
is changed, the peaks of the SPA shift whereby all curves share the static SPA as a limit. 
This results into the qualitative frequency dependence found by McNamara and Wiesenfeld 
for the delta function part of the spectrum of the double- well system (4). 

Roughly speaking, an improvement of the SPA under ferromagnetic coupling only occurs 
at d < 1, where the impact of the static SPA is not yet outperformed by the dynamical 
factor. This explains why there is no improvement at signal frequencies u >> 1: Since under 
ferromagnetic coupling 0(1 — 7) < 1 holds, d < 1 cannot be fulfilled in this case. At 00 = 10 
(Fig. 3) the SPA is already within line width of the limit p s /d 2 = a 2 {l — 7 2 )(u;T)~ 2 , where the 
impact of the dynamical factor prevails. This limit is basicly given by the modulation term 
in (15). Its non-monotonous temperature dependence is known to be no longer associated 
with a matching of time scales [10]. 

The static response, which is thus the origin of the coupling-induced enhancement of 
the SPA, can be interpreted as follows. Without coupling and for a given fixed signal 
the two-state elements prefer to be in the state of lowest energy (8) or highest barrier 
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(12), respectively. Since the signal is homogeneous, this state is the same for all elements. 
Without a signal but under ferromagnetic coupling neighbouring elements tend to be in 
the same state, too, although they do not favour a particular state +1 or —1. Taking now 
signal and coupling together, both effects add up. The tendency to find the elements in the 
state favoured by the signal grows and hence q s and the static SPA grow compared to the 
uncoupled case. 

With antiferromagnetic coupling, on the other hand, neighbouring elements prefer to be 
in opposite states. This counteracts the effect of the signal and leads to a decrease in the 
SPA. 

Finally, the decrease of the static SPA with growing temperature directly follows from 
(9). There the imbalance in the distribution of probability between the two states, which 
dependent on the view taken either differ in energy or barrier height, decreases as the system 
heats up: The mean response q s weakens due to increasing fluctuations. 

V. COUPLING AND SIGNAL-TO-NOISE RATIO 

The investigation of the SNR presented here basicly relies on numerical evaluations of 
(25). An explicit calculation of the power spectrum (21) already yields a rather complicated 
expression which does not lend itself to a detailed analytical study. 

At finite signal frequencies u> > the SNR displays the same qualitative dependence 
on temperature T and coupling parameter J as shown for the SPA in Fig. 2 and Fig.3, 
respectively. At low frequencies the SNR is enhanced under ferromagnetic coupling whereas 
at high frequencies it is not. Due to this close similarity we omitted the respective plots for 
the SNR. We note, however, that this similarity is by no means a trivial result, since the 
spectrum (21) is itself a non-monotonous function of T and J. 

The low-frequency SNR qualitatively reproduces the SNR behaviour found in a chain 
of next-neighbour-coupled overdamped double- well systems simulated in [8]. It is the first 
analytical confirmation of this behaviour. We note that this correspondence occurs although 
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the simulation in [8] was performed with strong forcing while the present model is studied 
within a weak-signal limit. As for the SPA, the high-frequency behaviour of the SNR has 
to our knowledge not been reported before. 

If the elements are not coupled, their SNR (3) does not depend on the signal frequency. 
This was already found in [10], if the signal-induced reduction of the continuous part of the 
spectrum is neglected, as (3) clearly shows. Under coupling this frequency independence of 
the SNR is lost. For ferromagnetic coupling numerical evaluations of (25) predict a decrease 
of the SNR with growing signal frequency. Only in the limit of low and high frequencies the 
SNR approaches constant values. In both cases the SNR formula (25) simplifies significantly. 
One finds 



where i?o is the simplified McNamara-Wiesenfeld SNR (3) of uncoupled elements. R s is the 
static SNR while Rhf represents the leading order term of its high-frequency expansion. 
Fig. 4 illustrates the typical qualitative behaviour of the SNR as a function of the signal 
frequency. At various frequencies it shows a set of SNR curves over temperature T and 
coupling parameter J, respectively. The dashed curves represent R s and Rhf, respectively. 

The apparent fact that the SNR decreases with growing signal frequency implies that its 
static value cannot be exceeded at any other signal frequency. Then it follows immediately 
from (32) that the coupling-induced enhancement of the SNR possesses an upper limit: It 
cannot be better than a factor four. This is a rigorous result for the general chain model and 
not limited to the prototype system, where a special choice of the parameters was made. Due 
to the simplicity of the chain model this result may well reflect a limit of more general nature 
for the improvement of a weak-signal SNR. Fig. 5 shows the static SNR over temperature 
T. It differs qualitatively from the static SPA which according to (24) has a monotonous 
behaviour. 




(32) 




(33) 
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Turning briefly to antiferromagnetic coupling one finds that the continuous part of the 
spectrum is insensitive to the sign of the coupling parameter J. Now even the static SNR 
decreases under coupling. As shown in Fig. 4 (bottom) it is found below the high-frequency 
SNR curve, i.e., here the SNR increases with growing signal frequency. Since at medium fre- 
quencies again a transition occurs between both curves similar to the transition in Fig.4(top), 
a second maximum in the SNR over T may emerge in this situation. 

The investigation of the impact of the coupling on the SNR can also be extended to 
the more natural situation, where the signal itself is embedded into noise. The question to 
address is, whether the improvement of the SNR, which was established for independent 
internal noise sources, can still be found with additional coherent external noise. (This 
problem does not occur for the SPA, where only the height of the signal peak is of interest.) 

To study this case it is assumed that the input spectrum consists of the previous signal 
peak described by ttHq5(H — oo) and a noise part N(Q). Within the weak-signal limit it was 
shown that the signal peak of the input spectrum leads to a respective peak in the output 
spectrum at signal frequency, only: No additional peaks at multiples of that frequency occur. 
Therefore, any additional contribution to the continuous part of the output spectrum at 
signal frequency can only arise from N(Q — uj). In analogy to (16) one finds that this 
additional contribution is given by (q/ H ) 2 N(u). 

Together with (25) the resulting SNR R noiS y reads, again in units of Hq, 

Rnmsy = HgS(u) + q*N(uy (34) 

It can be expressed in terms of the SNR R (cf.(25)) and the input SNR Ri nput = ix/N(oj) 
written in units of Hq, respectively. One finds 

1 1 1 



Rnoisy R Rinput 



(35) 



Rnoisy is thus a steadily growing function of the previously studied SNR R, whereby it cannot 
exceed the input SNR Ri npu t- (The latter is a meanwhile well-known result of linear response 
theory [18] which is in fact the limit we are taking here.) Since Rinput is constant at fixed 
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signal frequency u, a coupling-induced increase of R will lead to an increase in R no isy, too. 
The improvement of the SNR under coupling is thus preserved with external coherent noise. 
From (35) one can easily show that the maximum enhancement of R no isy compared to the 
uncoupled chain reaches the previously found factor four for vanishing input noise, only. In 
general, this factor is smaller approaching 1, if the input noise is so strong that the input 
SNR and hence Rnoisy go to zero. 

At the end we shall make a brief excursion to the collective response of the chain, which 
was mentioned in the introduction. We consider the following sum which involves the states 
of N elements: 

N 

M(f) = 5>(f). (36) 

i 

It is not difficult to show that the SPA of this new quantity is simply N 2 times the previ- 
ously studied SPA (23). For the SNR, however, this new situation is completely different. 
Previously the continuous part of the spectrum was determined from the auto-correlation 
function Ckk(t,r) of an element alone (cf. (18)). Now this spectrum will involve cross- 
correlation contributions Cjk(t,r) of different elements j ^ k, too. 

The new contributions to the continuous part of the spectrum change the qualitative 
behaviour of the SNR under coupling. We shall demonstrate this in the limit N — > oo. 
For this case Glauber calculated the unperturbed spectrum [11] which reads in its one-sided 
version 

s "<"> = ^£Sp - (37) 

Inserting this spectrum into (25), multiplied by N 2 due to the mentioned increase of the 
SPA, one finds for the new SNR per element 
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R M = -a^l- 1 2 5 2 , (38) 

which, of course, again reduces to the simplified McNamara-Wiesenfeld SNR (3) for vanish- 
ing coupling. 
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Clearly, the SNR is now a decreasing function of the coupling parameter and the sign 
of the latter is no longer important. Hence, there is a drastic difference in the SNR of 
single and collective response with respect to coupling. Of course, we have so far merely 
investigated the limits of (36), N — 1 and N — > oo, respectively. We expect that the SNR of 
the collective response of only a few elements still increases under coupling [16]. A detailed 
analysis of the present model with respect to these collective effects will be the subject of a 
further investigation. 

VI. SUMMARY 

In this paper we extended the two-state model of stochastic resonance introduced by 
McNamara and Wiesenfeld to a chain of infinitely many coupled two-state elements which 
are periodically modulated. The interaction of the elements was chosen in a way that 
the chain evolves according to Glauber's stochastic Ising model. In analogy to the work of 
McNamara and Wiesenfeld and based Glauber's results analytical expressions for the signal- 
to-noise ratio (SNR) and the spectral power amplification (SPA) have been obtained in the 
limit of weak modulations. Here both quantities refer to the response of a single element as 
part of the chain. 

Instead of approximating the dynamics of a particular coupled system by the chain model, 
we used the latter to build a prototype system which hopefully captures the essential features 
of an entire class of coupled stochastic resonators. To this end additional assumptions 
were made on the dependence of the model parameters on some noise intensity and signal 
amplitude, respectively. The prototype system was used to study the effect of the coupling 
on the response of a single resonator embedded into the chain. 

The results show that array-enhanced stochastic resonance occurs for ferromagnetic-type 
coupling in SPA and SNR. The qualitative features of the effects reproduce those previously 
found in coupled stochastic resonators. The simplicity of the chain model allowed for a 
detailed analytical investigation of the SPA. For the SNR the observed effects have been 
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confirmed analytically for the first time. In addition, it was found that the improvement 
of the single-resonator response compared to the response of the uncoupled resonator still 
occurs, if the signal is embedded into noise. For the SNR this improvement was shown to be 
limited by a factor four which is reached for vanishing input noise, only. A brief excursion 
into the collective response of N resonators, on the other hand, disclosed that coupling 
cannot improve the SNR, if N is very large. 

A closer look at the mechanisms behind the effects revealed that in the present model 
an improvement of the stochastic resonators under coupling is associated with the regime 
of quasistatic response. Since the model studied here possesses an upper bound to the time 
scale of its dynamics, the desired improvement is essentially restricted to sufficiently slow 
signals. For the SPA this improvement is based on a stronger tendency of the two-state 
elements to align in parallel, if signal and ferromagnetic coupling act together, compared to 
this tendency caused by the signal alone. The reason why the improvement of the SPA is 
lost, if the coupling is too strong, was found to be the slow-down of the system dynamics 
under coupling: It simply prevents the resonators from responding quasistatically. 

We would like to thank an anonymous referee whose valuable recommendations helped 
to shape the final version of this paper. 
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FIG. 1. Features of the unique maximum of the SPA p over temperature T and coupling 
strength J shown in dependence on the signal frequency u. Plotted are the position (T max , J max) 
of the maximum, its height p ma x, and the sharpness of the maximum expressed by the curvatures 
Ptt and pjj, respectively. In the uncoupled model (J = 0) the SPA has a maximum at To with a 
height pq. 
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FIG. 2. Typical qualitative behaviour of the SPA p at low signal frequencies for various values 
of coupling strength J (top) and temperature T (bottom). Dashed and dotted curves correspond 
to SPA limits discussed towards the end of this section. The dash-dotted curves represent the SPA 
at optimal coupling. 
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FIG. 3. Typical qualitative behaviour of the SPA p at high signal frequencies for various values 
of coupling strength J (top) and temperature T (bottom). 
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FIG. 4. The frequency-dependence of the SNR R. The upper dashed curves represent the static 
SNR. The lower dashed curves show the high-frequency expansion of the SNR. 
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FIG. 5. Static SNR R for various values of the coupling parameter J. The SNR at J = 10 is 
within line width of the amplification limit 4Rq- 
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